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Abstract 

The classical Stefan problem is reduced as the singular limit 
of phase-field equations. These equations are for temperature u and 
the phase-field ip, consists of a heat equation: 

ut + iipt = An, 

and a Ginzburg-Landau equation: 

eift = eAip - ^W'{ip) + l{ip)u, 

where f is a latent heat and fU is a double-well potential whose wells, 
of equal depth, correspond to the solid and liquid phases. 

When e —7- 0, the velocity of the moving boundary v in one dimension 
and that of the radius in the cylinder or sphere are shown as the fol¬ 
lowing Stefan problem, 


( Ut — An = 0 


^ 2 


du 

dn 


J r 


( u = -^[K-av]r 

where a is a positive parameter, [|^]r is the jump of the normal 
derivatives of u (from solid to liquid), and m = {2W{ip))^^‘^ dip. 

Since it is sufficient to describe the phase transition of single compo¬ 
nent by the phase-field equation, we analyze the phase-field equation, 

^ = o?Aip + f {ip), 

where a is a positive parameter, and f{ip) a function of the double¬ 
well potential and investigate whether the equation shows the Stefan 
problem or not. The velocity of the moving boundary in the cylinder 
and sphere are determined and the result of the simulation of the 
equation is also presented. Next, we consider the velocity of interface 
which depends on the temperature. The control of width of diffusion 
layer by the parameter of phase-field equations is investigated in order 
to realize the singular limit of phase field equations by the numerical 
method. 
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1 Introduction 


Phase-field equations are used in the field of phase transition. Recently, the 
relationships between the physical properties and the parameters of these 
equations have been determined and as a result, phase-held equations have 
been established as the position of the model for simulation of the phase 
transition behavior in engineering held. 

Phase-held equations are also applied to the double layer model of interface 
between huids for example, interface of liquid-gas and liquid-liquid contact. 
However in the cases of contact of gas-sold and liquid-solid, we should mod¬ 
ify these equations for the dihusion layer not to be double, but the width of 
dihusion layer tends to zero. 

The singular limit of the phase-held equations have been investigated the¬ 
oretically and numerically. Caginalp and a co-researcher [1-4] have studied 
the reduction of these equations in the singular limit of Stefan problems and 
Hele-Shaw problem. Furthermore H. Soner [11] showed that the solution of 
the phase-held equations converged to that of the Stefan problem. An asymp¬ 
totic behavior of solutions for the Stefan problem is calculated by Carslaw 
and Jaeger [5], Friedman [6]. In their study, the moving interface distance on 
the process of the solidihcation is a half-square of time, namely, this result is 
coincided with the dependence on time as the free moving interface length, 
with two points from the origin, of the phase-held equations. 

It is sufficient to describe the phase transition of single component by the 
phase-held equation. Therefore we analyze the equation and obtain the 
asymptotic behavior of the equation in order to investigate that the equation 
shows the Stefan problem. 

This paper consists as follows. Sec.2 is devoted to the review of studies by 
many researchers in order to make clear the relations among the theories con¬ 
cerning of this problem for applied mathematician, physicist and engineer. 
Sec.3 is a numerical example. 


2 Review on Mathematics of This Problem 

Caginalp and co-researchers reported that the phase-held equations con¬ 
verged to the equations for the Stefan problem as singular limit. These 
equations are for temperature u and the phase-held if, consists a heat equa¬ 
tion: 

ut + ift = Au, 
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and a Ginzburg-Landau equation: 

1 

e(pt = eA(f - 

€ 

where £ is a latent heat and hh is a double-well potential whose wells, of 
equal depth, correspond to the solid and liquid phases. 

When e —)■ 0, the velocity of the moving boundary v in one dimension 
and that of the radius in the cylinder or sphere are shown as the follow¬ 
ing Stefan problem and let r{t) be the interface separating the two regions 
Q(t) = {x : (p(t,x) = — 1 } or {99 = 1 }, 


Ut — Am = 0 , in G = {x : </:’(t,x) = ± 1 } 
1 du 

^ 2 (9n p 

7Tl 


where a is a positive parameter, [|^]r is the jump of the normal derivatives 
of u (from solid to liquid), and m = dip. On the other hand, 

lida [7] suggested that those equations did to the equations for the two phase 
Stefan problem by the same singular limit. These of the fundamental idea 
is based on the work about the model of “Caginalp type”. More precisely, 
these results can be derived by the method of matched asymptotic expansion 
between the inner expansion and the outer expansion. 

Soner reported that the solution of the phase-held equations converges to 
that of the Stefan problem, when e —)■ 0. Now, when e —)■ 0, the solution 
(nb if'') of 

r ui + = Au% 


converges to the solution (m, (p) of 


{ Mi - Am = -{h{p>))t 
u = ±^[i^-av]r 

where - 1 - stands for solidihcation and — melting, and i{(p) = h\p>) = a/ 2 W {p>). 
From these results, the moving velocity of interface is obtained from the limit 
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of the phase-field equations. The velocity of interface without the interface 
tension is 



t/ — —1— 

am 

when the temperature n = 0 on T, the velocity 

1 

V = ± — K 

a 

where a is a positive constant. 

3 Numerical Example 

We consider the following equations, 

(3.1) = €^Alp + f{ip). 

Suppose that the width of the interface between solid and liquid e is suffi¬ 
ciently small in comparison with the radius of solid R{t) as shown in Fig.l. 
Let us consider the time variation R{t). In the case that e/R 1, we can 
write 

(3.2) <I>(r,t) = <l>o(r ±i?(t)). 

Substituting (3.2) to (3.1), we get 

(3.3) ± m', = + ^^ 0 ) + /(<^o). 

where d is a dimensional number, 2 or 3. 

Since the equilibrium is assumed at the interface, the following equation is 
given, namely, 

(3.4) + 

holds. Here <hg is a large value at the neighborhood of r = i? and rapidly 
decrease to zero out of that neighbor. Hence we replace ‘ho/r to ^'q/R and 
as a result, from (3.3), we get 

(3.5) = 


{R{t)f = Rl + 2e^{d-l)t 


and 

(3.6) 
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Figure 1: R(t), Rq and x, ^/-coordinate 


where Rq is the initial radius. The relation of {Rq vs. t is good agreement 
with tendency of the exact solutions [5] of the Stefan problem. 

Next, we consider the velocity of interface which depends on the tem¬ 
perature. The control of width of diffusion layer by the parameter of phase- 
held equations is investigated in order to realize the singular limit of phase 
held equations by the numerical method. Karma and Rappel [8] in 1996 
derived parameters for the so-called thin-interface limit of phase-held equa¬ 
tions, where the interface thickness can be controlled to be small, and the 
classical interface conditions are satished for a hnite thickness. Their anal¬ 
ysis allowed for the hrst time fully resolved computations to be made for 
three dimensional dendrites with arbitrary interface kinetics, by Karma and 
Rappel [9] in 1998. 

We used their model to simulate the Stefan problem by the phase-held equa¬ 
tions. In this simulation, we gave the function as 


(3.7) 
and 

(3.8) 





du 

m 


Au + 


1 d(p 

2~m 


(3.9) -^ = A^ + W\^)-X{u) 

The right hand side of the equations is introduced to the parameter A. The 
parameter A is controlled by the thickness of dihusion layer between solid 
and liquid, and is similar to the one used by Caginalp and others (e.g. see 
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MacFadden, et. al. (1993) [10]). 

Now, we consider the steady state of one dimensional phase-field equation to 
compare the sharp interface model. The boundary conditions of the equations 
are 


(3.10) 


V = 


du 

dx' 


(3.11) 


Ui = -I3V, 


(3.12) u{+oo) = —6. 

where in the equation (3.11), subscript i stands for the interface. Using the 
equations, we obtained the solution of the steady state equations as 

(3.13) r = 

(3.14) u = exp[—Vx]—5 

in the liquid {x ^ 0) and n = 1 — 5 in the solid (x ^ 0). 

In one dimension, phase-held equations take the form 

(3.15) W = ll + F - ^’) + 


(3.16) 


du d'^u 1 dip 
dt dx"^ 2 dt 


The steady-state growth equations in the moving frame of interface, which 

+ p - + \u = 0, 


du ^ d‘^u V dp ^ 


yield 


(3.17) 

dx 

(3.18) 

V 


dx dx"^ 2 dx 


The solution of these equations with u subject to the far held boundary con¬ 
ditions determines the planer interface velocity as a function of undercooling. 
We numerically analyzed the equations (3.17), (3.18). Figure 2 and Figure 
3 show the prohles of p and U. The thickness of the interface is dependent 
on the parameter A. The thickness of p and U with A decreasing. This 
procedure of thin interface limit gives the solution of the Stefan problem as 
singular limit of phase-held equations. 
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4 Conclusion 


Rearranged the theorems of phase-field equations in singular limit of 
Stefan problem, we reviewed the previous work reported by mathematician. 
The velocity of the moving interface between solid and liquid was simply 
determined. 
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